Semiclassical propagator of the Wigner function 



O 

o 

(N 

(N 

> 

vn 
o 

oo : 
o - 

in 

o ■ 

Oh: 

c ■ 

=3 ; 



Thomas Dittrich,^ Carlos Viviescas,^ and Luis Sandoval^ 
^ Departamento de Fisica, Universidad Nacional, Bogota D.C., Colombia, 
Max Planck Institute for the Physics of Complex Systems, Nothnitzer Strafie 38, 01187 Dresden, Germany 

(Dated: February 1, 2008) 

Propagation of the Wigner function is studied on two levels of semiclassical propagation, one 
based on the van-Vleck propagator, the other on phase-space path integration. Leading quantum 
corrections to the classical Liouville propagator take the form of a time-dependent quantum spot. 
Its oscillatory structure depends on whether the underlying classical flow is elliptic or hyperbolic. 
It can be interpreted as the result of interference of a pair of classical trajectories, indicating how 
quantum coherences are to be propagated semiclassically in phase space. The phase-space path- 
integral approach allows for a finer resolution of the quantum spot in terms of Airy functions. 
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Quantum propagation in phase space has always been 
in the shadow of propagation in conventional (position, 
momentum) representations. Yet it is superior in various 
respects, particularly in the semiclassical realm: It avoids 
all problems owing to projection, such as singularities at 
caustics. Canonical invariance of all classical quantities 
involved is manifest. Boundary conditions are imposed 
consistently at a single (initial or final) time, thus remov- 
ing the so-called root-search problem and allowing for 
initial-value representations. Semiclassical approxima- 
tions to the quantum-mechanical propagator have pre- 
dominantly been seeked in the form of coherent-state 
path integrals Q,!^!^. Closely related are Heller's Gaus- 
sian wavepacket dynamics |^ and its numerous modifica- 
tions. By now, a broad choice of phase-space propagation 
schemes is available which score very well if compared to 
other semiclassical techniques. 

Almost all of these developments refer to the prop- 
agation of wavefunctions in some Hilbert space. Less 
attention has been paid to the propagation of Wigner 
and Husimi functions. They live in projective Hilbert 
space, i.e., represent the density operator and are bilin- 
ear in the wavefunction. Besides their popularity, they 
have a crucial virtue in common: An extension to non- 
unitary time evolution is immediate. This opens access 
to a host of applications that combine complex quan- 
tum dynamics, where a phase-space representation facil- 
itates the comparison to the corresponding classical mo- 
tion, with decoherence or dissipation: quantum optics 
and quantum chemistry, nanosystems in biophysics and 
electronics, quantum measurement and computation. 

By the scales involved, many of them call for semi- 
classical approximations. However, only few such stud- 
ies exist, for specific systems predominantly in quantum 
chaos [3, including dissipative systems 0, H- By con- 
trast, Ref. m discusses a new method, Wigner-function 
propagation analogous to the solution of classical Fokker- 
Planck equations. 

As a major challenge, any attempt to directly prop- 
agate Wigner functions requires an appropriate treat- 
ment of quantum coherences. As early as 1976, Heller 
argued that the "dangerous cross terms", i.e., the off- 



diagonal elements of the density matrix in the relevant 
representation, can give rise to a complete failure of semi- 
classical propagation of the Wigner function. Quantum 
coherences are reflected in the Wigner function as "sub- 
Planckian" oscillations They plague semiclassical 

approximations by their small scale and by propagating 
along paths that can deviate by any degree from classical 
trajectories. 

In this Letter, we point out how Heller's objections 
are resolved by considering pairs of trajectories as ba- 
sis of semiclassical approximations, and present corre- 
sponding expressions for the propagator of the Wigner 
function. The concept of trajectory pairs has been in- 
troduced in the present context by Rios and Ozorio de 
Almeida |Il[l3, 

albeit working in a strongly restricted 
space of semiclassical Wigner functions. We here give 
a general derivation of the propagator, independently of 
any initial or final states. 

Moreover, we go beyond the level of approximations 
based on stationary phase. Employing a phase-space 
path-integral technique, we construct an improved semi- 
classical Wigner propagator in terms of Airy functions. 
It resolves all singularities and contains the semiclassical 
approximations based on trajectory pairs as a limiting 
case. The interference patterns we obtain depend, up 
to scaling, only on the nature of the underlying classi- 
cal phase-space flow — elliptic vs. hyperbolic — and in this 
sense are universal. While living in projective Hilbert 
space, this result is superior to Gaussian wave-packet 
propagation in that it allows Gaussians to evolve into 
non-Gaussians. 

In order to fix units and notations, define the 
Wigner function corresponding to a density operator p 
as W{r) = J d^q' exp(-ip • q'/h) (q + q'/2| p |q - q72) 
where r — (p, q) is a vector in 2/-dimensional phase 
space. Its time evolution is generated by a Hamiltonian 
H{p, q) through the equation of motion {d/dt)W{r, t) = 
{H{r),W{r, t)}Moyai, involving the Weyl symbol H{r) of 
the Hamiltonian H. The Moyal bracket {.,.}Moyai 
converges to the Poisson bracket for ?i — > 0. As this 
equation of motion is linear, the evolution of the Wigner 
function over a finite time can be expressed as an inte- 
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gral kernel, W{r",t") = J d^^r' G{r" ,t" ir' ,t')W{r' ,t'), 
defining the Wigner propagator G{r" ,t";r' ,t'). For au- 
tonomous Hamiltonians, it induces a one-dimensional dy- 
namical group parameterized by i = t" — t' (in what fol- 
lows, we restrict ourselves to this case and use t as the 
only time argument). This implies, in particular, the 
initial condition G(r",r';0) = S{r" ~ r') and the com- 
position (Chapman-Kolmogorov) equation G{r",r',t) = 
Jd^frGir",r,t-t')G{r, r',t'). 

The Wigner propagator can be expressed in terms of 
the Weyl transform of the unitary time-evolution op- 
erator U{r,t) = /d-'^q'exp(-ip • qy/i) {q + q'/2\ U{t) 
q — q'/2), called Weyl propagator, as a convolution, 

G(r", r', t)= J d2-^'i?e^("""-^')^^C/*(f_, i)C/(f+, i), (1) 

with r± = (r' + r" ± R)/2. It serves as a suitable start- 
ing point for a semiclassical approximation, invoking an 
expression for the Weyl propagator semiclassically equiv- 
alent to the van-Vleck approximation [T3. ITHj |. 



cxp (iS'j(r, i)/?i — i^jTT/2) 
V|det(M,-(r,t7TT^ 



(2) 



The sum includes all classical trajectories j connecting 
phase-space points r^ , r" in time t such that r = f ^ = 
(r^ + Mj is the corresponding stability matrix, 

/ij its Maslov index. The action Sj {vj , t) = Aj (r^ , t) — 
H{rj, t) t, with Aj, the symplectic area enclosed between 
the trajectory and the straight line (chord) connecting 
and r" (hashed areas Aj± in Fig. 

Substituting Eq. (0) in Eq. leads to a sum over 
pairs j+, j-, of trajectories whose respective chord cen- 
ters fj± are separated by the integration variable R. 
Otherwise, the two trajectories are unrelated. A cou- 
pling between them, as expected on classical grounds, 
comes about only upon evaluating the R-integral by sta- 
tionary phase. Stationary points are given by r" — r' = 



(r"_ 



r^_l-)/2. Together with the conditions 



for the two chords, r' -|- r" ± R = r'j^ + r"_|_ , this implies 



(r;-+r;+)/2. 



)/2. (3) 



Stationary points are thus given by pairs of classical tra- 
jectories such that the initial (final) argument of the 
propagator is in the middle between their respective ini- 
tial (final) points (Fig. 13)). This does not require these 
trajectories to be identical! They do coincide as long as r' 
and r" are on the same classical trajectory, but bifurcate 
as r" moves off the classical trajectory rci(r',i) starting 
at r', if the dynamics is not harmonic. 

The resulting semiclassical approximation for the Wig- 
ner propagator is (dot indicating time derivative) 



G(r",r',t) ^ -jY: 



Af ^ 2cos(S'j(r",r',i)/n-/7r/2) 



y ^ ^|det(M,+ -M,_)| 
5,(r", r', t) = (f,+ - f,_) A (r" - r') + S,+ - 



(4) 



= / ds[f,{s) A R,(s) - i?,+ (r,+) + i/,-(r,_)] , (5) 
Jo 




FIG. 1: The reduced action (shaded) of the Wigner propaga- 
tor in van-Vleck approximation, Eq. Q , is the symplectic area 
enclosed between the two classical trajectories Vj±{t) and the 
two transverse vectors r'j^ — r'j_ and r"_|_ — r"_ (schematic 
drawing). The full central line is the classical trajectory 
rci(r',i), the broken line is the propagation path fj(r',t). 



with v,it) = (r,_(t) + rj+(t))/2, R,(t) = rj+(t) -r,_(0, 
and Sj± = Aj±{rj±,t) — Hj±{rj±,t)t. The reduced ac- 
tion Aj — Jq ds fj(s) A Rj(s) is the symplectic area en- 
closed between the two trajectory sections and the vec- 
tors r'^ — r'j_ and r"^_ — r"_ (Figs. 1112b ). In general, 
Eq. 0}, as a function of r", describes a distribution that 
extends from the classical trajectory into the surround- 
ing phase space, forming a "quantum spot" (Fig.Oa) with 
a characteristic oscillatory pattern that results from the 
interference of the contributing classical trajectories. In 
general, it fills only a sector with opening angle < 27r 
(Fig. Et), where the sum contains two trajectory pairs 
(four stationary points). Outside this "illuminated area" , 
stationarity cannot be fulfilled, that is, the "shadow re- 
gion" is not accessible even for mean paths rj{t). The 
border is formed by phase-space caustics along which 
there is exactly one solution (two stationary points). As 
r" approaches the classical trajectory rci(r',t) starting 
at r', from the illuminated sector, the two solutions j — , 
= j+ coalesce so that Afj_ — > Mj+, and Eq. becomes 
singular. If the potential is purely harmonic, all mean 
paths coincide with rci(r',t), and the classical Liouville 
propagator, G(r",r',t) = (5(r" — r"j(r', t)), is retained. In 
all other cases, Eq. though based on the van-Vleck 
propagator, reflects the structure of stationary points of 
the action including third-order terms, with one pair of 
extrema and one pair of saddle points. It is formulated in 
terms of canonically invariant quantities related to clas- 
sical trajectories and thus generalizes immediately to an 
arbitrary number of degrees of freedom. The propaga- 
tion of Wigner functions defined semiclassically in terms 
of Lagrangian manifolds ^| is contained in Eq. Q as a 
special case. 

We are now able to resolve Heller's objections 0: If the 
two trajectories j — , j+ are sufficiently separated and the 
potential is sufficiently nonlinear, then (i), the propaga- 
tion path r{r' ,t) can differ arbitrarily from rci(r', t), and 
(ii), the phase factor in exhibits sub-Planckian oscil- 
lations. They would couple resonantly to corresponding 
features in the initial Wigner function, generating a sim- 
ilar pattern in the final Wigner function around the end- 
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FIG. 2: Classical building blocks entering the Wigner propa- 
gator according to Eqs. I)4I5|I . for a stable trajectory starting 
at r' = (0.636, 0) near the minimum of the cubic potential 
V(g) = 0.329q'^ — 0.69g (panel a), (b) Classical trajectory 
rci(r',i;) (black line), a pair of auxiliary trajectories Vj±{i) 
(green lines) and corresponding propagation path f j (r' , t) (red 
dashed line). The yellow target pattern is the grid of aux- 
iliary initial points around r', parameterized by polar 
coordinates. Propagated classically over time t, it deforms 
into the turquoise pattern around r". The red/blue cone is 
formed by midpoints f" = (r"_ + r"^)/2 that correspond to 
extrema/saddles of the action. Its boundaries form caustics 
separating the region accessed by two midpoints r" from the 
unaccessible rest, (c) Enlargement of the area around r", 
indicating the number of trajectory pairs that access each re- 
gion. 



point of the non-classical propagation path. In this way, 
quantum coherences are faithfully propagated within a 
semiclassical approach. 

Equations I|4I5|I translate into a straightforward algo- 
rithm for the numerical calculation of the propagator 
(Fig. 121): (i) Define a local grid (e.g., in polar coordinates) 
around the initial argument r' of the propagator, identi- 
fying pairs of auxiliary initial points r^^, r^_|_ with r' in 
their middle, (ii) Propagate trajectory pairs rj±(i) clas- 
sically, keeping track of the symplectic area Aj between 
them, (iii) Find the amplitude and phase contributed by 
each trajectory pair and associate them to the final mid- 
points f". They constitute a deformed cone, projected 
onto phase space (Fig. 13;). Its "lower" ("upper") surface 
(red (blue) in Fig. |3;) corresponds to pairs of extrema 
(saddles) of the action, respectively: (iv) Superpose the 
contributions of the two surfaces, after smoothing ampli- 



tude and phase over midpoints r" within each of them. 

The caustics in Eq. result from applying stationary- 
phase integration in a situation where pairs of stationary 
points can come arbitrarily close to one another. Since 
the underlying van-Vlcck propagator admits only up to 
quadratic terms in the phase, we seek a superior ap- 
proach, corresponding to a uniform approximation. It 
is available in the form of a path-integral representation 
of the Wigner propagator 16], in close analogy to the 
Feynman path integral. 



G(r",r',t) ^l^jDrJ Di^e-'^»"■>'{^»/^ 



Two paths, r{t) and R(<), have to be integrated over. 
The former is subject to boundary conditions r(0) = r' 
and r{t) = r", the latter is free. The path action is 

5({r},{R})= / ds[r(s)AR(s) 

+Hw{r{s) + R(s)/2) - Hvj{r{s) ~ R(s)/2)] . (7) 

Equation is recovered upon evaluating the path- 
integral representation in stationary-phase approxima- 
tion: Defining r±{t) = r{t) ± R(t)/2, with boundary 
conditions analogous to Eq. (PJ, and requiring station- 
arity leads to the Hamilton equation of motion for r± (t) : 
We again find pairs of classical trajectories that straddle 
the propagation path as stationary solutions. 

We will now include cubic terms in the action, with 
respect to variations of the path variables. To keep 
technicalities at a minimum, we restrict ourselves from 
now on to a single degree of freedom and to Hamiltoni- 
ans of the standard form H(p,q) = T{p) -\- V{q), where 
T(p) = /2m while the potential V{q) may contain non- 
linearities of arbitrary order. With this form, Hw{r) 
coincides with the Hamiltonian function "quantized" by 
merely replacing operators with classical variables. As 
the path integral readily allows to treat time-dependent 
potentials, chaotic classical motion remains within reach. 

Expanding the action |7J around r{t) — rci(r',t) and 
R(t) = {P{t),Q{t)) ~ 0, there remain only linear terms 
in P and linear and cubic terms in Q. Evaluating the R- 
sector of the path integral thus results in an Airy spread- 
ing of the propagator, with a rate ^ V"'{qci{t)), in the 
p-direction. It is superposed to the classical phase-space 
flow around the trajectory, i.e., rotation (shear) if it is 
elliptic (hyperbolic). As a consequence, a spot of the full 
phase-space dimension develops. Scaling p = (r?,^) = 
{^'^/'^p^^i^^/^q), with ^ = T"{pc\)/V"{qc\), we express 
the linearized classical motion as a dimensionless map. 



M{m) 



cos (j){t) 
sin 0(i) 



— sin 4>{t) 
cos (j){t) 



(8) 



These maps form a group parameterized by the angle 
(t>{t) = J* ds VT"(p'=i(s))V"(g=i(s)). It is real (imagi- 
nary) if the linearized dynamics is elliptic (hyperbolic). 

This allows to evaluate also the r-sector of the path 
integral. Transforming the Wigner function to Fourier 
phase space, W^-j) = {FW){'y) = (27r)~^ / d^r exp(-i7A 
r) W{r), and the propagator accordingly, G = FGF^^, 
we obtain (7' = (a',/3')) 



G(7",7',t)='5(7"-M(0")7') 
exp(-i(— 



'3 + a2ia'^f3' + ai2a'/3'2 + ^/3'3) ) . (9) 



(6) 



The coefficients ajk = /g dstT(s)(sin0(s))-'(cos(/)(s))'^ de- 
pend on where along the classical trajectory how much 
quantum spreading - a{t) = {^{t)f/'^h^V"'{qcx{t))/S is 
picked up and thus on the specific system and initial con- 
ditions. The Fourier transform from Eq. back to the 
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FIG. 3: Quantum spot replacing the classical delta function 
on a stable (elliptic) trajectory near the minimum of a cubic 
potential as sho-wn in Fig.|5Ji, at f = 1.8 {4> ~ 27r/3). Panel (a) 
shows the exact quantum result for the Wigner propagator, 
(b) and (c) are semiclassical approximations based on Eqs. @ 
and respectively, all for h — 0.01. Frames coincide with 
that of Fig.|2J:. (d) Quantum spot, according to Eq. Q, for an 
unstable classical trajectory near the maximum of the same 
potential, at t = 1.0. Crosses mark the classical trajectory. 
Colour code ranges from red (negative) to blue (positive). 



original Wigner propagator can be done analytically, af- 
ter transforming the third-order polynomial in the phase 
to a normal form ■ 

The internal structure and the time evolution of the 
quantum spot described by Eq. (|5J) are qualitatively dif- 
ferent for elliptic and hyperbolic classical trajectories 
(real and imaginary (p, resp.). In the elliptic case, the 
spot is a periodic function of (p. In particular, it collapses 
approximately to a point whenever = 2Z7r, I integer. 
Close to these nodes, it shrinks and grows again along a 
straight line in the p-direction, reflecting the fact that for 
short time, the quantum Airy spreading ~ t^^^ outweighs 
the classical rotation ~ t. Only sufficiently far from the 
nodes, while rotating around the trajectory by (j)(t)/2, 
the one-d. distribution fans out into a two-d. interference 
pattern formed as the overlap of the bright (oscillatory) 
sides of two Airy functions, with a sharp maximum on 
the classical trajectory (Fig.|3|D). In comparison with the 



corresponding exact quantum- mechanical result (Fig.|2K) 
for the quantum spot, obtained by expanding the prop- 
agator in energy eigenstates |0| , the path- integral solu- 
tion resolves the caustics far better than Eq. Q (Fig.|3;). 
The hyperbolic case is obtained replacing trigonometric 
by the corresponding hyperbolic functions. As a result, 
along unstable trajectories there are no periodic recur- 
rences as in the elliptic case; the spot continues expand- 
ing in the unstable and contracting in the stable direction 
(Fig. 131). Isolated unstable periodic orbits embedded in 
a chaotic region of phase space exhibit a degeneracy of 
the Weyl propagator It allows to account for scar- 
ring in terms of the Wigner propagator |0| . 

We have obtained a consistent picture of incipient 
quantum effects in the Wigner propagator, both in the 
van-Vleck approach and in the path-integral formalism: 
(i) for anharmonic potentials, the delta function on the 
classical trajectory is replaced by a quantum spot extend- 
ing into phase space, (ii) its structure shows a marked 
time dependence, qualitatively different for elliptic and 
hyperbolic dynamics, (iii) it exhibits interference fringes 
arising as a product of Airy functions, (iv) it can be ex- 
pressed in terms of canonically invariant quantities as- 
sociated to pairs of underlying classical trajectories, (v) 
within each level of semiclassical approximation used, the 
propagator retains its dynamical-group structure. Open 
issues include: extension to higher dimensions and to 
higher-order terms in the action, performance in the pres- 
ence of tunneling, application to unstable periodic or- 
bits and implications for scars, trace formulae, and spec- 
tral statistics, regularization of the ballistic nonlinear a- 
model, semiclassical propagation of entanglement, and 
generalization to non-unitary time evolution. 
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